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Abstract Conditions in the protosolar nebula have left their mark in the com¬ 
position of cometary volatiles, thought to be some of the most pristine material 
in the solar system. Cometary compositions represent the end point of processing 
that began in the parent molecular cloud core and continued through the collapse 
of that core to form the protosun and the solar nebula, and finally during the 
evolntion of the solar nebnla itself as the cometary bodies were accreting. Disen¬ 
tangling the effects of the various epochs on the final composition of a comet is 
complicated. But comets are not the only source of information about the solar 
nebula. Protostellar disks around young stars similar to the protosnn provide a 
way of investigating the evolution of disks similar to the solar nebula while they 
are in the process of evolving to form their own solar systems. In this way we 
can learn about the physical and chemical conditions under which comets formed, 
and about the types of dynamical processing that shaped the solar system we see 
today. 

This paper summarizes some recent contributions to onr understanding of both 
cometary volatiles and the composition, structure and evolution of protostellar 
disks. 

Keywords protostellar disks • solar nebula • comets • chemistry 
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1 Introduction 


Cometary volatiles are some of the most pristine of solar system materials, having 
remained relatively unprocessed since the comets formed in the inner solar system 
some 4.2 billion years ago. Consequently they retain in their composition signatures 
of the chemical and physical conditions under which they formed. Observations of 
cometary comae provide the mixing ratios or relative abundances of various species 
and have demonstrated that these can vary greatly between comets, e.g.,|A’Hearn| 


et al. (20121 and Mumma and Charnley (2011) and references therein, with little 


correlation with dynamical family. The origin of these compositional variations is 
not well understood but must lie in the physical and chemical processes that were 
active at the time and location at which the cometary volatiles were formed. To 
understand these processes requires not only a study of the comets themselves but 
also of protostellar disks similar to the one from which our solar system formed. 
Observations of protostellar disks allow for the important disk processes to be 
investigated while they are occurring and can test models of the evolution of the 
early solar system. 

Some key questions that can potentially be answered by interdisciplinary in¬ 
vestigations involving both protostellar disk modeling and observations, and the 
observations of comets are: 


1. Can synergy between protoplanetary disk modeling and compositional studies 
of comets inform us about the chemical history and processes in the early solar 
disk? 

2. What is the origin of the diversity of observed cometary compositions? 

3. Where did comets form in the solar nebula? 

4. Are D/H ratios in volatile molecules and spin temperatures cosmogonic indi¬ 
cators for comets? If so, what can they tell us about the formation conditions 
in the early solar system? 

5. What is the range of D/H ratios in cometary water and can we infer the degree 
to which comets contributed volatiles to the early terrestrial planets? 

6. Can cometary volatile abundances be used to constrain the extent of radial 
and vertical mixing in the early protoplanetary disk? 

This paper summarizes some recent developments in our understanding of 
cometary compositions and their relationship to the conditions under which the 
comets formed. The paper is arranged as follows. Section covers some obser¬ 
vational and modeling studies of protostellar disks. Section discusses fraction¬ 
ation processes in protostellar disks including using deuterium fractionation as a 
tracer of the links between comets, interstellar chemistry and nebular chemistry. 
Section considers the origin of the nitrogen deficiency in comets and Section 
discusses spin temperature as a cosmogenic indicator in comets. Section [^presents 
a possible explanation for the observations of neutral sodium cometary tails. Fi¬ 
nally, Section [T] summarizes this paper and briefly discusses the potential impact 
of Rosetta on our understanding of the formation conditions of comets. 


2 Overview of the structure and composition of protostellar disks 

Protostellar disks form during the star formation process. A molecular cloud core 
collapses under gravity forming a protostar in the center. As material continues to 
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fall inwards conservation of angular momentum results in the formation of a disk 
around the central protostar. In the early stages, the disks are gas-dominated with 
masses ~ 1% of the mass of the protostar. They rotate about the protostar with 
a Keplerian velocity prohle. Disks contain a wide range of physical conditions, 
ranging from cold (~ 10 - 20 K in the outer disk midplane) to hot (with gas 
temperatures up to several thousand K in the surface layers). Irradiation conditions 
also vary, with the surface layers experiencing the full force of the stellar UV field, 
but with little radiation reaching the midplane. Consequently there are a wide 
variety of chemical conditions depending on location within the disk. The star and 
disk formation process, and the chemistry and dynamics of disks are discussed in 
more detail in recent reviews by Dutrey et al. (2014); Pontoppidan et al. (2014) 
and Turner et al. (2014). 

Disks can be very turbulent places. Turbulence is difficult to detect directly but 
evidence for its effects can be seen in the presence of small dust grains high above 
the disk midplane. Without turbulence these grains would quickly be removed from 
these high altitudes by settling (Dullemond and Dominik 2004). The dynamical 
motions also drive collisions between grains, which can lead to both grain growth, 
eventually forming planetesimals, and to creation of new small grains or dust. 
With grain growth, the grains will eventually reach a size at which their motions 
decouple from those of the gas and they begin to settle under the influence of 
gravity to the midplane, where they can undergo further collisions and growth to 
form planetesimals. The removal of grains from the surface layers allows the stellar 
UV to penetrate deeper into the disk, changing its composition. Turbulence can 
also drive chemical changes by bringing together material processed in different 
parts of the disk and conversely chemistry can affect the dynamics by controlling 
the ionization levels. Consequently dynamics and chemistry strongly influence each 
other (see Section 2.3.2). 

In recent years protostellar disks have been an area of active research, with their 
molecular content, dust properties, and physical conditions being studied in detail. 
The observational work has been coupled with simulations that can interpret the 
observations and make predictions about those regions of the disks which are 
not amenable to direct observational study. This synergy of the observations and 
modeling work has revealed a wealth of information about the chemistry and 
physics of disks, and in turn has informed our understanding of the early history 
of our own solar system. This section reviews some of the recent observational and 
modeling work related to protostellar disks, concentrating on the early stages of 
evolution when the disk is still gas-dominated (the T Tauri phase). 


2.1 The structure of a protostellar disk 

Protostellar disks are characterized by strong radial and vertical gradients in tem¬ 
perature and density (Figure[^. The stellar UV field controls the temperature and 
hence the chemistry throughout much of the disk. The penetration of the UV into 
the disk is governed by the opacity which is mainly provided by dust absorption. 
In the surface layers, where the UV field is unattenuated, the photoelectric effect 
can efficiently heat the gas to much higher temperatures than the dust. Close to 
the midplane, where the UV field is reduced, the dust and grain temperatures are 
coupled. 
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Fig. 1 The physical structure of a typical T Tauri disk. Shown are the density (left), gas 
temperature (center) and dust temperature (right). The density and dust temperatures are 
calculated by a hydrodynamical model of Paola D’Alessio with a similar mass and surface 
density profile to the minimum mass solar nebula. The gas temperature is determined via 
thermal balance and is coupled to the dust temperature throughout most of the disk, but is 
much hotter in the surface layers due to the photoelectric effect. The black contour line on the 
density plot indicates the likely location of the dead zone, determined by the location where 
the ionization level x(e) < 10“^^. In this region the turbulent motions driven by MRI will be 
severely inhibited. 


These gradients in physical conditions lead to a wide variety of chemical envi¬ 
ronments (Figure]^. Models have found that outside of the water ice line the disk 
can be divided vertically into three chemically distinct layers. At the surface, the 
disk is dominated by the effects of UV. Ices are desorbed (either thermally or by 
photodesorption) and efficient photodissociation means that the gas is mainly com¬ 
posed of atoms and their ions. In the midplane most molecules are frozen out onto 
the surfaces of dust grains. Between these two layers is a region where the grains 
are warm enough to allow at least some molecules to thermally desorb and which 
is shielded enough that photodissociation is inefficient, allowing the molecules to 
survive once in the gas. It is this warm molecular layer that is detected in most 


gas phase observations of the outer regions of protostellar disks (see Section 2.2). 

As the distance to the protostar decreases the disk warms up. In the midplane, 
the ice mantles will begin to desorb, with different molecules being released at 
different radii, depending on their binding energies - the most volatile species being 
released at larger radii than the less volatile ones. This results in a series of ice lines. 


which mark the transition from ice to vapor for each molecule (see Section 2.41. The 


icelines are curved, and in the outer disk are almost horizontal with the molecule 


existing as ice below this line and as gas above it (see Pontoppidan et ah] 2014). 
Once inside the water iceline (within a few AU of the protostar) the molecules are 
all in the gas phase. 

Figure also illustrates some of the main physical processes that are active in 
disks. In particular, turbulence is a crucial process, enabling dust coagulation and 
driving chemical changes by mixing material processed under different physical 
conditions within the disk. What drives this turbulence is still a matter of debate 
but one of the main candidate processes is the magneto-rotational instability (MRI; 


e.g., Balbus and Hawley|1991 Hawley and Balbus|I991 Hawley et al.|1996 ) which 
can efficiently transport angular momentum in ionized accretion disks. A minimum 
level of ionization is required (in the min imum mass solar nebula the ionization 
fraction, x(e), must be > 10“^^ at 1 AU; Inutsuka and Sano 2005). This can be 


difficult to achieve in the shielded midplane regions of a protostellar disk, where 
the high column densities prevent UV and cosmic rays from penetrating. Hence 
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Fig. 2 The physical and che mical environment of a 1 - 5 M yr old protostellar disk around a 
Sun-like star. Reprinted from |Henning and Semenov| l |^13[ l with permission. @(2013) Amer¬ 
ican Chemical Society. 


there are regions in the disk where the MRI is not active. These are known as 
“dead zones” and they are predicted to occnr throughont much of the midplane in 
the planet-forming region. The approximate location of the dead zone in the disk 
shown in Figure is indicated by the black contour line in the density plot. 

The dependence of MRI-driven turbulence on the ionization means that there 
is a close link between the chemistry and the dynamics. The chemistry determines 
the ionization level and therefore the strength of the turbulence, but the dynamical 
motions themselves can drive chemical changes through mixing. This in turn affects 


the ionization and consequently the dynamics (e.g., 

InutsTika and Sano|2005 

Ilgner 

and Nelsen|2008 

Ilgner and Nelson|2006 

Turner et al.|2007 

1. Hence the evolution 


of chemistry and dynamics are inextricable linked. 


2.2 Observations of protostellar disks 


Protostellar disks have been observed at wavelengths from the UV to the mil¬ 
limeter. Each wavelength can probe a different region of the disk. The majority of 
information about the molecular composition of disks comes from the infrared and 
the (sub)-millimeter. The disk is mostly optically thick in the infrared, therefore 
these observations are restricted to the warm gaseous region inside of the iceline 
or to the surface layers. At these wavelengths gaseous molecules can be observed 
in both absorption and emission. Ices can also be observed in absorption if illumi¬ 
nated by the central star in an edge-on disk, e.g., Pontoppidan et al. (20051, but 


the spectra can be difficult to interpret due to contributions from material outside 
of the disk. Near- and mid-infrared observations are most sensitive to the inner 
regions of disks, while the far-infrared can probe the surface layers at larger radii 
(see Henning and Semenov (2013). 
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Spitzer, operating in the mid-infrared, has shown that inside of the iceline 
the disk is chemically rich, with the terrestrial planet-forming region containing 
organics such as C 2 H 2 , HCN, as well as CO, CO 2 , water and OH (e.g., Carr and 


Najita 2008 

Salyk et al. 

2008 

Pontoppidan et al. 2010 


the result of a rich gas phase molecular chemistry, in a region where temperatures 
range from a few hundred kelvin up to 1000 K, and densities are high (> 10® 


cm“®). 

Spitzer observations of the gaseous inner disk have also been used to make 
inferences about the conditions beyond the iceline. For example, [Najita et al.| 
(2013) found the HCN/H 2 O ratio inside of the iceline to be correlated with disk 
mass, suggesting that this ratio may reflect changes to the C/0 ratio induced by 
the formation of icy planetesimals outside of the iceline. Outside of the iceline 
much of the oxygen is tied up in water ice, whereas the carbon is mainly in the gas 
phase in the form of more volatile molecules such as CO. Planetesimals outside of 
the iceline will form from grains coated with water ice. Since the planetesimals do 
not move about the disk as easily the small grains and the gas, this will trap the 
oxygen outside of the iceline and the material that moves inwards to the inside 
of the iceline will be enriched in carbon. Consequently the gaseous inner disk will 
be enriched in organics compared to water. Since more massive disks are likely to 
form planetesimals more efficiently than their lower mass counterparts they will 
have a higher C/0 ratio inside of the iceline. 

In the far-infrared, Herschel has detected a number of molecules including cold 
water in the outer disks of TW Hya (Hogerheijde et al. |2011[ ) and DC Tau ( |Podio| 


et al.|20l3). In both cases the water emission comes from a region of the disk where 


grain temperatures are low enough for this molecule to be expected to be frozen 
onto dust grains. Its presence in the gas can be explained by photodesorption 
driven by the stellar UV field (Willacy and Langer 2000 Dominik et al.j 2005l. 


Hogerheijde et al. (2011) found that the predicted flux from photodesorbed water 


to be higher than observed, and suggested that the answer to this discrepancy could 
be grain coagulation. As the grains grow they decouple from the gas motions and 
sink towards the midplane, taking their ices with them. This removes water ice 
from the molecular layer, leaving less to be photodesorbed, and thus reducing the 
gas phase abundance. 

Emission lines from some molecules in the disk can be optically thin in the 
(sub)-millimeter and so in theory these wavelengths can probe the whole of the 
disk. In practice, however, resolution constraints restrict these observations to > 
30 - 100 AU from the protostar, although ALMA allows observations to within a 
few AU of the closest protostars e.g. (ALMA Partnership et al. 2015 ). Observations 
by both interferometers and single dish facilities have detected the emission from 
molecules in the warm molecular region and their distributions can be mapped 
using interferometers. . Relatively simple molecules have been observed including 
CO (and its isotopologues), H 2 CO, CS, C 2 H, HCN, CN, H NC, DCN, HC 3 N, c - 
C 3 H 2 , N 2 H+, HCO'*", DCO+, H 2 D+ (for a recent review see Dutrey et al.||2014 |. 
The fractional abundances relative to H 2 are generally lower than observed in 
molecular clouds because of the combined effect of depletion in the cold dense 
midplane and photodissociation in the surface layers. Both processes reduce the 
gas phase molecular abundances. 

In addition to the composition of the disk, the detection of molecules can also 
provide information about the physical conditions under which they exist and the 
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chemical processes that are active (Table[^. Observations of molecular deuteration 
are especially important since they trace the thermal history of the disk (see also 
Section 3.11. Deuterated molecules, with their lower column densities, also provide 
a tracer of the gas in the midplane, e.g., ALMA observations of DCO^ have been 


used to trace the location of the CO iceline in a Herbig AeBe disk (Mathews et 


al. 


2013) and H 2 D+ is a potential tracer of the ionization level in the midplane 


( Ceccarelli et al.]|2004 ). 


Process 

Molecules 

Temperature 

CO, i^co, Ci® 0 , H2CO, H2O, OH, H2 

Density 

CO, H2CO, HCO+, HC3N 

Ionization 

HCO+, N2H+, CH+ 

Velocity 

CO, CS, HC3N 

Photodissociation 

CN, HCN, HNC, OH, H2O, C2H 

Photodesorption 

H2O, OH 

Grain surface chemistry 

H2CO, complex organics, CO2, C2H2 

High temperature chemistry 

Complex organics 

Deuteration 

HD, DCO+, DON, H2D+ 


Table 1 Molecules observed in protostellar disks and the processes they trace. 


2.3 Protostellar disk modeling 


Chemical models of protostellar disks vary in their complexity, but they all predict 
the three-layer structure beyond the iceline as shown in Figure]^ The chemical 
simulations require as their basis, a background disk model to provide the den¬ 
sity, temperature and UV held as a function of location in the disk. These can 
be provided by hydrodynamical models of a 1-| -1D steady-state Q-dis k model in 

digosl p' 


vertical hydrostatic equilibrium such as those of 


D’Alessio et al. 


(see also 


Figure]^. These physical disk models are computationally intensive and are are 
generally run independently of the chemical simulations. Current chemical mod¬ 
els are based on detailed chemical kinetics and can contain hundreds of species 
(in both the gas and the ices) and thousands of reactions. These models include 
the ionizing effects of UV photons, cosmic rays. X-rays and short-lived radionu¬ 
clides such as ^®A1. Some models have also focused on isotopic chemistry such as 


deuterium (e.g., 

Aikawa and Herbst||1999 

Woods and Willacy||2007 

Willacy and 

Woods||2009 Kavelaars et al. 

2011 Furuya et al.|2013 

Albertsson et al.||2014 

1 or 

carbon and oxygen isotopes ( 

Woods and Willacy|2009 

see also Section 

3.2 

1 . The 


simplest models assume that the conditions in the disk did not change with time 


and ignore dynamical motions, although a few also consider advection e.g., Willacy 


et al. 

19981; 

Aikawa et al. 

(1999) 

Woods and Willacy 

(20071 

2009 

); Willacy and 

Woodf 

(2009 

). More recent models have included processes such as grain growth 


http;//www.cfa.harvard.edu/youngstars/dalessio/ 
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and settling, e.g., Aikawa and Nomura (2006); Fogel et al. (2011 

; Vasyunin et al. 

(2011 

; Akimkin et al. (2013 

), dynamical motions, (e.g. Ilgner et al.|2004 |Willacy 

et al. 

2006t Semenov et 

al. 

2006[ Semenov and Wiebe||2011[ |Ilgner and Nelson 

2006 

Turner et al.||2006 Aikawa||2007 Hersant et al.|2009 Heinzeller et al.|2011 1 

and hydrodynamical evolution of the disk structure itself (e.g. I 

lodson-Robinson 


2.3.1 Formation of complex organics in protostellar disks 


Of particular interest to the chemistry of comets is the work of Walsh et al. (20141 


who modeled the possible formation of complex organics throngh ice chemistry in 
protostellar disks. Figurej^shows the range of abundances they find in ices located 
outside of 20 AU in a solar-type disk. The large molecules form by sequential 
addition of atoms to existing mantle molecules, e.g., ethylene glycol can form by 
adding hydrogen, oxygen and carbon atoms to CO: 

CO ^ HCO HC 2 O OCCHO 
OCCHO ^ CHOCHO ^ HOCH 2 CHO ^ (HOCH2)2 


(Charnley and Rodgers 2008[ see Figure]^. The models of Walsh et al. suggest 
that chemistry in disks will greatly increase the abundance of these large molecules 
relative to what can be created in the parent cloud. The predicted abundances 


relative to H 2 O are consistent in many cases with observations of comets (Bockelee- 
Morvan et al. 2004 Crovisier||2006 Crovisier et al. 2004). If these molecules can 
be desorbed they will also be present in the disk gas, and should be detectable by 
ALMA providing a test of the model. 
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Fig. 3 Complex organics can be formed by chemistry in the ice mantles of dust grains in 
a protostellar disk. This figure shows the predicted range of abundances of some of these 
molecules r elative to water ice (red lines) compared to those derived from com etary comae 
(blue lines; [Bockelee-Morvan et al. |2004 |Crovisier]|2006| |Crovisier et al.||2004| |. The green 
asterisks are the results from the initial niolecular cloud rnodel used as the input for the disk 
model. From jWalsh et aL] j |2014[ |. 
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2.3.2 The effects of dynamics on the chemistry 


There are strong links between the disk chemistry and its dynamics. The tempera¬ 
ture, density and radiation field in a disk control the chemistry, but the chemistry 
(and grain evolution) in turn controls the dynamics of the disk through the ioniza¬ 
tion structure. Models have found that better agreement with observations of the 
molecular layer in the outer disk can be obtained if vertical mixing is included, e.g.. 


Willacy et al. 

(2006) 

Semenov et al. 

(2006), and more 

sophisticated models have 

now considered 2-D mixing (e.g.. 

Heinzeller et al.|2011 

Semenov and Wiebe|2011 


Albertsson et al. 20141. Mixing smooths out chemical abundance gradients but also 


brings together species that might not otherwise exist in the same location, e.g., 
in the inner disk, vertical mixing combines OH and H 2 in a region where the tem¬ 
perature is high enough for the neutral-neutral reaction, OH -|- H 2 = H 2 O -|- H to 
occur ( Heinzeller et al.]|2011 |. As a result the abundance of H 2 O increases. 

In addition to affecting molecular abundances, mixing also changes isotopic 
ratios. One problem with the non-mixing models is the high deuteration levels pre¬ 
dicted in the comet formation zone, e.g., Willacy and Woods (2009). These reflect 


the high levels set in the parent molecular cloud, whereas those detected in comets 
tend to be somewhat lower. As D/H ratios track the thermal history of the material 
this suggests that at least some of the cometary ices formed at warmer temper¬ 


atures than are found in the interstellar medium (ISM), e.g., Meier and Owen 
( 1999[ ). Figure [^compares the observed molecular deuteration in cometary water 
with the predictions of several models. The non-mixing models uniformly over¬ 
predict HDO/H 2 O, but the addition of mixing (either radial or vertical) brings 
the ratio into closer agreement with the observations ( Bonev et al.||20i4 ). 

With vertical mixing this reduction in the deuteration of water is achieved 
by the destruction and reformation of ice mantles. When dust grains are mixed 
vertically they can reach the surface layers where high temperatures and/or UV 
fields desorb their ices. The water molecules are dissociated and the resulting 
oxygen, hydrogen and deuterium atoms travel back down towards the midplane 
into more shielded regions where they reform molecules. The molecular deuteration 
depends on the formation temperature. Since the disk is warmer than the 10 K 
parent molecular cloud core, the reformed water ice has a lower D/H ratio. Hence 
vertical mixing can reduce the initial interstellar D/H ratio i n water ice ~ 10~^ ) 
to the level detected in comets (~ 10“^) (Figu rean d see also Furuya et al. 20131. 
Willacy and Turner (2015) and Furuya et al. (2013) also predict radial variations 


in the deuteration of water. A radial gradient in HDO/H 2 O is also found in the 
(radial) mixing models of Kavelaars et al. (20111. 


2.3.3 Grain evolution 


Grain-grain collisions in protostellar disks lead to efficient grain growth and ulti¬ 
mately result in the formation of planetesimals. As the grains grow they will de¬ 
couple from the gas motions and sink towards the midplane. This removes grains 
from the upper layers of the disk, reducing the opacity and allowing UV to pene¬ 
trate deeper, which increases photodissociation and heating of the newly irradiated 
layers, and moves the molecular layer closer to the midplane. An additional effect 
of grain growth is a reduction in the freezeout rate, which depends on the surface 
area of the grains. Figureshows the effect on the chemistry of CO in the models 
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Water D/H measurements in comets 
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Fig. 4 Current disk models can only acco unt for the observed D/H ratios in comets if they 
include mixing. The non-mixing models of|Willacy and Woods| (|2009|l and T. Millar (private 
communication) predict D/H ratios that are consistently higher tha.n the observations. |A1-| 


_ on ) 

|bertsson et al.| ||2011|l treat the evolution of the molecular cloud to the hot cor e stage only and 
predict /^)water > 10“^. Turbulent mixing is included in the disk models of |Kavelaars et ah] 
|20TTJ| (radial mixing) and IWillacy and Turner] ||2015|| (vertical mixing at 5 AU) and these 
produce T)/Yi water iu agreement with t he observations. Furuya et al.| j2013|| p redic t similar 
D/H r atios to |Willacy and Turner] l|2015|l . Reference s: (1) Eberhardt et al. Jl|1995[|, (2)|Hartogh | 


et al.| ||2011[ ), (3) I Bockelee-Morvan et al.l j2012||, (4)|Lis et a" 1 |2013|, (5)|Bockelee -Morvan et 
al.| p99 8|l (6a-d ) |Meier et al.H|1998|l and|(Jrovisier et al.| | |2004 i, (7 )|Hi 


,|l|1995|l, 

|Bockel6 


^ utsemekers et al.|l|2008|| 

|Biver et al.| 1 |2006[ |) (9) |Viilanueva et' al.| 1 |2009[ |. Adapted from [Bonev et al.| ( |2014| |. 
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Fig. 5 The CO abundance distribution in a protostellar disk, showing the effects of grain 
growth. The highest abundances are shown in red. Grain growth reduces the opacity in the 
surface layers, allowing UV to penetrate further and pushing the molecular layer closer to the 
midplane. With 1 micron dust (left) the molecular layer is much thicker, and reaches closer 
to the midplane than in the model with 0.1 micron dust (center). On the right are the results 
from an evolving dust model where the grains are allowed to both grow and sediment. In this 
case the molecu lar layer is thicker tha n in the small grain model, and it moves closer to the 
midplane. From [Vasyunin et al.| l [2011[ |. (c)AAS. Reproduced with permission. 
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of Vasyunin et al. (20111 (see also Aikawa and Nomura||200'6 Walsh et al.||2014| . 
Vasyunin et al. found that CO, CN and SO, as well as the ratio of C 2 H 2 /HCN, 
should be good tracers of the grain growth. 


2.4 Gas-dust chemistry and icelines in protoplanetary disks 


Icelines in the early nebula are dehned as radii at which specihc molecules undergo 
a shift in abundance from the gas phase to the condensed phase. There are three 
of specific interest related to H 2 O, CO 2 and CO and these occur at radii where 


the temperature is approximately 150 K, 47 K and 20 K respectively (Oberg et 


al. 20111. Recent observations by ALMA have allowed icelines to be identified 
in protoplanetary disks. For example, in HD 163926 | Qi et aL||2011 de Gregorio- 
Monsalvo et al.|2013 Mathews et al. 20131. The icelines associated with H 2 O, CO 2 , 
and CO are regions of dramatic chemical change and these could be reflected in 
the molecular composition of the gas and ice that is transported and mixed, both 
vertically and radially, in the evolving disk. 

The chemical signature of such an environment may be imprinted in cometary 
and meteoritic composition if the Sun formed in a stellar cluster containing at 


least one massive star (Adams 


2010 


Mumma and Charnley 20111. Nearby OB 


stars strongly affect the positions of the CO 2 and CO icelines and the chemical 
structure of the disk ( Walsh et al. 12013 ). The CO iceline occurs at tens of AU from 
the central protostar and is most accessible to observations at millimeter and sub¬ 
millimeter wavelengths. The physical conditions here most closely resemble those 
in the densest molecular cloud cores and so the related gas-dust chemistry could 
also be similar ( Henning and Semenov|2013 ). This chemistry can lead to formation 
of many complex organic molecules on cold dust grains, as well as potentially large 
isotopic fractionation in both gas phase and grain-surface reactions. 

For example, outside the CO iceline, CO-rich ices on dust grains can undergo 


(tunneling) addition reactions with H (and D) atoms accreted from the gas (Tielens 
1983[ Charnle^|1997 |. Atoms of O, C and N will also accrete and these processes 
can lead to rich organic chemistry (Figure[^ |Charnley|1997 Charnley and Rodgers 
2008 Herbst and van Dishoec^|2009[ ). Many of the proposed hydrogenation se¬ 
quences have now been demonstrated in surface chemistry experiments (Watanabe 
and Kouchi|2008 Theule et al.|2013 ) and such an ice chemistry explains the pres¬ 
ence of the distinctive and inter-related suite of organic molecules detected around 
massive and low-mass protostars (in hot cores and hot corinos, respectively), where 
the ice mantles have been evaporated into the hot gas ( [Herbst and van Dishoeck 
2009). Figure 1^ indicates that many of the relevant organic molecules could also 
be present in comets and that, if correct, ketene, ethanol, propynal, acrolein, pro- 
pionaldehyde and propanol should be searched for in suitably bright apparitions. 
Several of the organic compounds formed in the scheme of Figure [^ may be de¬ 
tectable in protoplanetary disks with ALMA. Vertical transport of ice-covered dust 
to the less well-shielded layers of the disk, as well as (inward) radial transport, can 
lead to these organic compounds being released by evaporation and/or photodes¬ 
orption. This is probably the origin of the H 2 CO detected in HD 163296, which 
lies outside of the CO iceline at about 160 AU ( Qi et aL]|2013b ). 
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Fig. 6 Grain surface reaction scheme forming organics via atom addition reactions to CO on 
cold dust. Molecules in red are detected both in comets and in interstellar/protostellar sources, 
whereas those in blue are found only in the latter. Species in black are not yet detected in 
space. The addition of deuterium atoms rather than hy drogen atoms to grain species can result 
in high levels of molecular deuteration. Adapted from|Charnley and Rodgers||^08||. 


2.5 Two distinct reservoirs of water ice in the protosolar nebula 


Formation scenarios of the protosolar nebula (PSN) invoke two main reservoirs 
of ices present in the disk midplane and that took part in the production of icy 
planetesimals (see Figure [^. The first reservoir, located in the inner region of the 
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PSN, contains ices (dominated by H 2 O, CO, CO 2 , CH 4 , N 2 and NH 3 ) originat¬ 
ing from the ISM, which, due to their near solar vicinity, were initially vaporized. 
The ice vaporization distance never exceeds ~ 30 AU from the Sun, depending on 


the total source luminosity and characteristics of the collapsing cloud (Chick and 
Cassen|[l997). With time, the decrease of temperature and pressure conditions al¬ 


lowed the water in this reservoir to condense at ~150 K (at typical nebula pressure 
conditions) in the form of microscopic crystalline ice (Kouchi et al. 19941, leav¬ 
ing negligible water in the vapor phase to condense at lower temperatures where 
amorphous ice would be expected. It is postulated that a substantial fraction of 
the volatile species were trapped as clathrates during this condensation phase as 


long as free water ice was available within 30 AU in the outer solar nebula (Mousis 


et al. 20001 and there was enough time to overcome the kinetics of clathration. 


On the other hand, the remaining volatiles that were not enclathrated (due to 
the lack of available water ice or a low kinetics of clathration) probably formed 


pure condensates at lower temperatures in this part of the nebula ( Mousis et al. 
2012a|b). The other reservoir, located at larger heliocentric distances, is composed 


of ices originating from the ISM that did not vaporize when entering into the disk. 
In this reservoir, water ice was essentially in the amorphous form and the other 


volatiles remained trapped in the amorphous matrix (Owen et al. 
and Bar-Nun|[2005 ). 


1999 


Notesco 


Interestingly, it has been shown recently that amorphous ice could form from 
photodesorption and freeze-out of water molecules near the surface layers of the 
PSN (Ciesla 2014). In these conditions, the transport of the icy grains through¬ 
out the outer solar nebula would lead crystalline ice to be lost and reformed as 
amorphous ice. In turn, amorphous ice generated on the disk’s surface could be 


transformed to crystalline as it migrates through the protoplanetary disk (Ciesla 


2014). This cycling scenario needs to be further investigated if one wants to assess 


its role in the formation or destruction of the two original icy reservoirs. Also, the 
different transport mechanisms must be simultaneously taken into account at the 
different stages of the disk’s evolution to depict the abundance of water ice in the 
two reservoirs as a function of the heliocentric distance. For example, as described 
in Sec. 2.3.2, it has been shown that vertical mixing can decrease the abundance 
of water inside the ~15 AU of the disk. In contrast, it has also been proposed 
that the outward vapor diffusion induced by the presence of the iceline with its 
subsequent condensation in a narrow location (with a width around 0.5 AU) would 
deplete the region inside the iceline of vapor down to subsolar values, and increase 
substantially the ice abundance at the iceline location by factors larger than 10 
times the protosolar value ( Ali-Dib et al.|[^014 ). In addition, recent works depict¬ 
ing the chemical evolution of both solids and gases from the pre-stellar core to 


the protostar and circumstellar disk phases (Visser et al. 


2009 


2011 


Harsono et 


al. 2013) should be taken into account to update the representation of the two 


water ice reservoirs present in the PSN. In these cases, midplane ices may have 
experienced multiple desorption events during the evolution of the disk. Never¬ 
theless, if one follows this classical picture, icy solids that formed at heliocentric 
distances less than ~30 AU mainly agglomerated from a mixture of clathrates and 
pure condensates, whose ratio depends on the amount of available crystalline wa¬ 
ter and its clathration efficiency ( Mousis et al.|2009b ). In contrast, solids produced 
at higher heliocentric distances (i.e., in the cold outer part of the solar nebula) 
were formed from primordial amorphous ice originating from the ISM. Therefore, 
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Fig. 7 Two reservoirs of ices in the PSN. A first reservoir (cold zone), located at distances 
higher than ~30 AU, is constituted from amorphous ices coming from ISM. The second reser¬ 
voir (hot zone), located within the ~30 AU, is made from volatiles initially in the form of 
amorphous ices that were transported from ISM towards the inner and hot part of the disk. 
When reaching regions with temperatures higher than ~150 K, these ices vaporized. During 
the cooling of the disk, volatiles located in the inner 30 AU condensed again but in crystalline 
forms, including both pure crystalline ices and clathrates (see text). 


depending where the outer solar system bodies formed, clathrates may have been 


agglomerated in comets (Marboeuf et al. 

2010 

201 l| 2012), and in the building 

blocks of the giant planets (Gautier et al. 

[2001 

Alibert et al.||2005a|b Mousis et 

al.|2009b 2012a) and in their surrounding satellite systems (jLunine and Stevenson 

1985 Mousis|2004 Mousis and Gautier|2004 Mousis et al.|2009a). Regardless the 


possible presence of clathrates, there is today indirect evidence that comets formed 
from crystalline water ice because their deuterium-to-hydrogen ratios measured in 
H 2 O are substantially lower (~6-12 times the protosolar value) than the most 
deuterium-rich primitive meteorite (~35 times the protosolar value), whose level 
of deuteration is expected to be close to the one found in the ice infalling from 
the presolar cloud (Kavelaars et al. 2011 Hartogh et al. 2011[ 


Ceccarelli et al. 


2014). This difference can only been explained by the isotopic exchange between 


the deuterium-rich water and molecular hydrogen in the gas phase (Kavelaars 
et al. ]SI [Hartogh et al.|[2011[ |Ceccarelli et al.||2014| ). Hence, in order to enable 
this reaction, amorphous ice entering the disk from ISM must have been vaporized 
prior recondensing again in crystalline form, otherwise the isotopic exchange did 


not occur. 


3 Fractionation in protostellar disks 

3.1 The deuterium fractionation: The link between comets and the earliest 
phases of the solar system and interstellar chemistry 


The origin of the volatile ices found in comets is a central issue for understanding 
the formation and early evolution of the solar system. There are three possibilities. 
First, one view is that all the ices are comprised of pristine interstellar molecules 
(e.g. Irvine et al.|1996 ). There is indeed a strong similarity with the many known 
interstellar molecules believed to form on interstellar grains, although there are 


differences in the relative abundances of the volatile ices (Mumma and Charn 


ley 2011). Until recently this idea was widely considered to be untenable (see 
and Cleeves et al. 2014a). Second, at the other extreme, it has been 


3.1.2 


Section 

advocated that any interstellar molecules experienced total obliteration and ‘chem- 
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ical reset’ in the disk (Pontoppidan & Blevins [2014 1. In this case all the volatile 
molecules observed in comets must have been formed by nebular chemistry. 

Comets and meteorites contain both high-temperature processed materials, 
such as crystalline silicates, and matter that apparently formed at very low temper¬ 
atures, as indicated by the isotopic enrichments measured in D and (Mumma 
& Charnley 2011). These characteristics are most naturally explained by the large- 
scale turbulent transport of dust and gas from different nebular chemical environ¬ 
ments into the comet-formation zones ( Ciesla and Charnley|2006 ). The recognition 
that non-equilibrium models of nebular chemistry are essential, particularly in the 
comet-forming regions, has motivated many detailed studies of protoplanetary 
disk chemistry, as recently reviewed by Henning and Semenov (20131 (see also 
Section 2.3.2). Turbulence leads to an outward radial diffusion of material from 
the hot inner nebula and an inwards advection of material from the cooler outer 
nebula. It also leads to vertical mixing and thus to midplane material experiencing 
the enhanced UV and X-ray fluxes of the upper disk layers (Glassgold et al. 1997). 
Hence, in this third partially-mixed scenario some specific volatile molecules and 
isotopic ratios could be remnant interstellar matter, or derive from warm and cold 
regions of the nebula (as reviewed by Pontoppidan et al. ]Ml4l ). Deciding between 
such molecular origins has been problematic. 


3.1.1 The molecular deuteration process 


Molecular deuteration provides a potential clue to the history of cometary volatiles. 
Although deuterium in the Universe is a tiny fraction with respect to the hydrogen 
(D/H=[2.535±0.05]X10“^: Pettini & Cooke 2012), its abundance in interstellar 
molecules has an enormous diagnostic power, because it may be much larger that 
the elemental D/H ratio. This generally occurs in cold (< 30-50 K) environments 
and in trace species, namely molecules with abundances smaller than ~ 10“®. 

The first key point to keep in mind is that, in cold gas, the major deuterium reservoir 
is HD and only a small fraction of deuterium is locked in other species. An exception 
to this rule is represented by the water, whose abundance can reach ~ 10““^ and, 
consequently, water might trap a substantial fraction of deuterium. 

The second key point is that the deuterium atoms are “extracted” from the HD 
molecules because the cosmic rays that hit H and H 2 create H^ ions, which in turn 
react with HD and form H 2 D'^ ions. 


HJ -h HD 


H 2 D+ + B.2 + 220K 


( 1 ) 


(see also Figure]^. The reverse reaction has a small energy barrier, which cannot 
be overcome at low (< 20-30 K) temperatures, so that the H 2 D+/HJ abundance 
ratio becomes larger than the D/H elemental ratio. The kinetics of Reaction 
at low temperatures is also strongly influenced by the spin state of the reacting 
hydrogen molecules, i.e., whether they are predominantly in the ortho or para 
states. The higher internal energy of ortho-}i 2 means that the reverse process 
can be driven efficiently if the H 2 ortho/para ratio is large. Thus, ort/io-H 2 acts 
as a “poison” for interstellar deuteration which requires that most of the H 2 be 
converted to para-Tl 2 (see Pagani et al.|20lT I. Reaction is the first step towards 
the deuteration of the trace species, as H 2 D"'“ reacts with all neutral molecules 
and passes, once every three times, the D atom to the product species. This is 












The composition of the protosolar disk and the formation conditions for comets 


17 



20141. 


Fig. 8 Schema of the molecular deuteration process in cold gas i Ceccarelli et al. 
Molecular deuterium occurs through three basic steps: 1) formation of ions from the 
interaction of cosmic rays with H and H 2 ; 2) formation of H 2 D"^ (HDJ and DJ) from the 
reaction of (H 2 D"*' and HDJ) with HD; 3) formation of other D-bearing molecules from 
reactions with H 2 D'^ and DJ) in the gas phase (Step 3a) and on the grain mantles 

(Step 3b). 


schematically shown in Figure A similar process can also occur through the 
CH 3 + ions at slightly higher temperatures. 

In practice, therefore, a molecular deuteration enhanced with respect to the 
D/H elemental abundance is a thermometer of the gas temperature. 

The third and final key point to keep in mind is that neutral molecules can freeze-out 
into the grain mantles and, in this way, they can be transmitted through the different 
stages of the star and planet formation. A sketch of this overall process is shown in 
Figure Therefore, ice species formed during the very first phases might be con¬ 
served and “reappear” in the gas at much later stages: for example when comets 
and planets are formed. Conversely, the molecular deuteration in comets and plan¬ 
ets can help us to understand in what epoch the relevant species was formed, at 
least in principle. In this sense, molecular deuteration can be thought of as an 
Ariadne’s thread which links the various phases of solar system evolution. For a 
more detailed review see Ceccarelli et al. |2014|. 


3.1.2 Deuteration in the solar nebula 

Deuteration in disks is driven mainly via the reaction of H 2 D"'“ ions (Reaction[^. 
At low temperatures the reverse process is inhibited and this leads to enhanced 
deuteration in gaseous molecules through ion-molecule isotope exchange reactions. 
Elevated atomic D/H ratios also pertain and lead to signihcant deuteration in 
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FROM A DIFFUSE CLOUD TO A SUN + PLANETARY SYSTEM 
THE DEUTERIUM FRACTIONATION HISTORY 

1- PRE-STELLAR CORE: cold (<10 K) & dense (=10=cm 3) gas 
FORMATION AND CONDENSATION OF HIGHLY DEUTERATED 
MOLECULES 


2- PROTOSTAR; warm (>100 K) & dense gas 

SUBLIMATION OF HIGHLY DEUTERATED MOLECULES 


3- PROTOPLANETARY DISK: 
cold, warm, & dense gas 
FORMATION AND CONDENSATION 
OF DEUTERATED MOLECULES 




4- PLANETESIMAL FORMATION : grains agglomeration 
& migration 

STORAGE OF DEUTERATED ICES AND REPROCESSING 


5- PLANET FORMATION AND COMETS/ASTEROIDES RAIN ON EARTH; 
planet migration, small bodies scattering; terrestrial oceans formation 
DELIVERY OF DEUTERATED MOLECULES 



Fig. 9 Schematic summary of the different phases of the formation of a solar-typ e star and 
planets, highlighting the formation and fate of deuterated molecules (adapted from|Ceccarelli| 
|et al.| ( (m4] l. Likely, the same process gave birth to the solar system. 


grain-surface reactions, such as those depicted in Figure and in the surface 
reduction of accreted O atoms to form HDO. Condensation of CO leads to en¬ 
hanced abundances of gaseous H 2 D+ and to s ignificant deuterium f ractionation 
beyond the CO iceline, as traced by DCO^ (Mathews et al. 2013). Enhanced 
D/H ratios in simple molecules, such as water, ammonia and methane, can be 
produced both in ion-molecule reactions and grain-surface additions of H and D 
to accreted atoms (Ceccarelli et al. 2014). At the CO iceline the grain-surface 
chemistry depicted in Figure 10 can lead to very high D/H ratios and multiply 
deuterated molecules ( Charnley et al.||1997 Ceccarelli et al. 2014). Observations 
show that simple deuterated species like DCN and by DCO"*" are present in disks 
(Qi et al. 2008) and we would expect that more complex fractionated organics, 
e.g., HDCO, D 2 CO, CH 3 OD and CH 2 DOH, to also be present. Between the CO 
and CO 2 icelines (see Section 2.4), large deuterium enrichments can persist in 
organic molecules since ion-molecule deuteration involving C 2 HD"'" and CH 2 D+ 
can persist up to higher temperatures (~ 70 K and ~ 50 K) than H 2 D+ (~ 30K) 
( Millar et al.|198'9 ). It is possible that the D/H ratios of several organic molecules 
can become re-fractionated even in the warm inner nebula, prior to the condensa¬ 
tion of acetylene and methane, both of which will occur outside the CO 2 iceline. 
This means that deuteration can still occur in cold gas that is mixed radially 
inwards, and be initiated in warmer gas coming from the inner nebula (Willacy 
2007 Albertsson et al.]|2014 |. 


However, the assumption of interstellar rates of cosmic-ray driven ionization 
(^) in the cold regions of the outer midplane has recently been called into question. 


Cleeves et al. (2013) showed that the effects of either deflection by stellar winds 


or magnetic field topology (through the magnetic mirror effect) could significantly 
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Fig. 10 Grain surface deuteration pathways associated with methanol formation from hydro¬ 
genation of CO ice. The probabilities for each reaction pathway, P^, assume a statistical steady 
state and are related to the atomi c D/H ratio of the ac creting gas and the energy barrier for 
addition of an H or D atom (from [Charnley et al.|1997^. 


attenuate the penetration of cosmic rays into magnetized protoplanetary disks, 
and substantially reduce ( by several orders of magnitude at the midplane. In this 
picture the only sources of ionization close to the midplane are scattered X-rays 
and the decay of short-lived radionuclides. In fact, Dolginov and Stepinski (19941 
had previously pointed out that inhomogenous random magnetic fields could also 
significantly inhibit the transport of cosmic rays throughout the disk and lead to 
very low values of C- 


Cleeves et al. (2014a I developed this scenario further by pointing out that effec¬ 


tively extinguishing cosmic ray ionization from the nebular midplane would greatly 
reduce the efficacy of Reaction and consequently the generation of enhanced 
molecular D/H ratios by nebular ion-molecule chemistry. They demonstrated that 
if the enhanced interstellar HDO/H 2 O ratios initially available were lost, then neb¬ 
ular chemistry could not regenerate those currently measured in the solar system, 
concluding that as much as 50% of the Earth’s oceans and perhaps all cometary 
water is of interstellar origin. Low values of C will also reduce the atomic D/H ratio 
in the disk midplane and so reduce the D/H ratios and the level of multiple deuter¬ 
ation in the molecules formed in the scheme of Figure[^ This may account for the 


upper limits on HDCO, CH 2 DOH and CH 3 DOD in comet Hale-Bopp ( Crovisier et 
ah]|2004 1 relative to the large observed interstellar D/H ratios measured in these 


molecules, which could have been lost or modified upon incorporation into the 
nebula. Although these molecules could be partially reformed in the cold outer 
nebula, this process, and any related deuteration, may be less efficient than in 
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molecular clouds. Cometary CH 3 OH/ H 2 O ratios are typically ~ 1 — 4%, whereas 
5 — 15% are typical for the ISM, and can be as high as 40% (Mumma and 


Charnley 2011). Cleeves et al. (|2014b) also modeled the general protoplanetary 


disk chemistry for a variety of assumed ( values and identihed molecular ions that 
may eventually allow the ionization structure to be determined. If these spatial 
distributions are observable with ALMA they could provide valuable information 
on the actual role of cosmic rays in nebular chemistry. 


S.i.S The message from molecular deuteration 

Figure shows the molecular deuteration in the different objects of the solar 
system and in objects which are believed to eventually form a solar-like system. It 
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Fig. 11 D/H ratio in molecules observed in solar system objects and in objects which are 
believed to eventually form a solar-like system. The filled symbols refer to water measures 
whereas open symbols refer to organic matter. The different objects are marked (note: CC 
means Carbonaceous Chondrites, lOM Insoluble Organic M atter, and SOM Solubl e Organic 
Matter). References to the various measures can be found in|Ceccarelli et al.|(|2014^. 


gives US the following major messages: 


— All reported solar system objects, including Earth, possess a molecular deuter¬ 
ation larger, by more than a factor 10, than the D/H elemental ratio of the 
environment where the Sun was born, and that is represented by the line 
marked with “Proto Solar Nebula”. This implies that the process of molecular 
deuteration was at work and it was at least partially inherited from the earlier 
phases of the solar system formation era. 
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There is a general trend, with the objects that formed earlier being also the 
ones with the largest molecular deuteration, in agreement with the scenario 
that the process started in very cold conditions. This implies that there was 
only a partial transmission between the subsequent steps of Figure Besides, 
earlier objects in the plot also correspond to material whose distance from the 
central heating source is larger. Therefore, the trend also tells us that there was 
not a general and substantial remixing of the different regions (outer and inner, 
where molecular deuteration is different because of the different temperatures). 
Exactly how much material was passed through each step or remixed remains, 
however, unanswered. 

The organic matter is systematically more enriched in deuterium atoms than 
the water. This is true in protostellar objects, comets and meteorites. However, 
while this result is robust for the protostellar objects and meteorites, where sev¬ 
eral measurements of the deuteration have been made, it is less certain for the 
comets, where so far the only organic species for which molecular deuteration 
has been measured is HCN. The reason for this dichotomy in the molecular 
deuteration is completely unclear. In protostellar objects, it is believed that 
this reflects the different phases of formation of water and organics, with water 
being formed very early when the gas is not very cold, while the organics were 
synthesized later from frozen CO. It is possible that a similar process is at 
the origin of the different deuteration in meteorites and comets. Alternatively, 
they may conserve a “core” of the protostellar molecular deuteration. 
Molecular deuteration in a given molecule has a strong temperature dependence 
at low temperatures. This is true for both gas phase and icy grain surface for¬ 
mation. In ice mantles the higher binding energy of deuteration atoms at < 
20 K leads to significant enhancements of deuterated species that form via hy¬ 
drogenation reactions (see also Figure [T^. These low temperatures correspond 
to what modelers predict for the cold midplane in the giant planet formation 
zone of the protoplanetary disk. Measurements of the molecular deuteration in 
comets provide a potential tracer of the thermal conditions under which they 
formed. For further discussion see Bockelee-Morvan et al. (this volume) 
Molecular deuteration in a given molecule has a strong temperature dependence 
at low temperatures. This is true for both gas phase and icy grain surface 
formation. In ice mantles the higher binding energy of deuteration atoms at 
< 20 K leads to significant enhancements of deuterated species that form via 
hydrogenation reactions (see also Figure 101. Measurements of the molecular 
deuteration in comets provide a potential tracer of the thermal conditions 
under which they formed. For further discussion see Bockelee-Morvan et al. 
(this volume). 


3 . 1.4 A model of the deuteration of water in the early stages of the evolution of the 
solar nebula 

Molecular deuteration potentially provides a means of tracing the history of solar 
system material and its links to the parent molecular cloud. In this section we 
present the results of a specific model of the evolution of the deuteration of water 
during the collapse of the parent molecular cloud, since the composition of the 
protosolar nebula depends in part on the chemical history of this core (Figure [T^. 
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The evolution of the core has been modeled recently by Aikawa et al. (20121; 
Wakelam et al. (2014); Taquet et al. (2014) and Tornow et al. ( 2014a|b ). 


Element abundances, temperature and density of inter-core matter 


Model of the quasi-stationary cloud core 


Chemical abundances, radial temperature and density distribution of core 


Model of the collapsing cloud core 


Chemical abundances, temperature and density distribution in disk range 

/•••» 

<> 

Fig. 12 Flowchart illustrating the evolution of the two precursor stages (the quasi-stationary 
cloud core and the collapsing core) to the solar nebula. Evolution in these stages will affect 
the composition of material incorporated into the solar nebula and hence the composition of 
planetary bodies that form within it. 


In these models the cloud core is initially described as a Bonner-Ebert sphere 
(Aikawa et al. 20081 which is embedded in a constant density gas-shell wiht a 
visual extinction of up to 4 magnitudes ([Wakelam et al.||2014), and irradiated by 


an extern al UV field. It is allowed to evolve for between 1 and 6 Myrs ( Tornow et 
al.|2014a|. Both the gas phas^and grair]^ chemistry is followed. Since the core is 


only partially shielded, photo-chemistry has a clear effect on the radial profile of 
the relative water abundance in the outer core region (Figure [l3|. At the end of 
the quasi-stationary core stage the temperature of the gas and dust decreases to 
8 - 9 K in the center for boundary values between 13 - 15 K. Consequently, water 
ice builds up on the cold grains, reaching to a fairly constant relative abundance 
of ~ 10“"^ within 10® - 10® years. In contrast. Figure 13 shows a clear radial 
gradient for the gas phase water abundance. Relatively low ratios (about 0.01 - 10 
ppm) of gas to dust phase abundance were obtained for the outer core range with 
up to 25 ppm in the inner region. Caselli et al. (2012) determined higher ratios 
due to a higher efficiency of water photodesorption; a consequence of the larger 
grain size and smaller photodesorption yield used in the Tornow et al. models. If 
the deuterium chemistry is included according to the method of [Albertsson et al.j 
(2013) the temporal and radial distribution of the gas phase HDO/H 2 O ratio 
can be derived. It reaches a maxiumum of about 0.02 in the medium shielded 
core region, and a minimum of 10”^ in the low shielded outer boundary region 
(Figure 14). In contrast, the dust phase ratio of HDO/H 2 O is clearly larger and 
only varies between O.OI and O.I. Its maximum is located in the well shielded inner 
core region at very early times (~ (2 -3) x 10® years). 


^ Nahoon_public_aprll.tar.gz at http://kida.obs.u-bordeauxl.fr/mod6ls/ 
^ disk_chemistry_OSU08ggs_UV.zip at http://www.mpia.de/homes/semenov 
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MRN Rad = l 


Zubko Rad = l 



MRN Rad = 3 


Zubko Rad = 3 


l.Oe-Oa 


l.Oe-09 


l.Oe-10 


l.Oe-11 


l.Oe-12 



Fig. 13 Relative gas phase abundance of ater versus time and radius (1 kAU = 10^ AU). The 
intensity of the standard interstellar radiati on field (IRSF) is sc aled by Rad = 1, 2, 3. Two 
grain distrib utions are used. MRN refers to |Mathis et al.| (|1977| and Zubko refers to |Zubko| 
|et al.|(|2004|. The mean grain size is 5.5 X 10“^ cm and 3.5 X 10“® cm for MRN and Zubko 
respectively. The cosmic ray ionization rate is tak en to be 10~^^ s~^. The related gas and 
grain chemistry is modeled using a combined code jTornow et al.|2014a[). 


The collapse of the cloud core forms a protostar surrounded by an optically 


thick disk and an outer envelope | Saigo et al. 


2008 


Schonke and Tscharnuter 


20111 . 


During this stage the chemical evolution can be followed along a Lagrangian path 
from the outer core boundary to the center of the core (Tornow et al. 2014bI. 
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Figure [T^ shows the HDO/H 2 O ratio along each member of the Lagrangian path 
set. Water formed in the quasi-stationary core stage is transported into the region 
of terrestrial planet formation. Signatures of large-scale, supersonic water inflows 
have been seen with Herschel iMottram et al. 20131 who determined that these 
flows occur on the core scale (< 10"^ AU), e.g., IRAS 15398 and L1527 or at least 
on the envelope scale (> 3 x 10^ AU), e.g., NGC1333-IRAS4A. Based on their 
data Mottram et al. suggest an outside-in collapse, in agreement with the model 
shown in Figure [Tsl 



l.OOOe+00 

l.OOOe-01 

l.OOOe-02 

l.OOOe-03 

l.OOOe-04 

l.OOOe-05 


2000 4000 6000 

R [AU] 


8000 


10000 


Fig. 14 Temporal and radial HDO/H 2 O distribution i n th e gas phase (Tornow et al. 2014b I 
for the standard ISRF and MRN grain size (see Figure [13^. 


The gas phase deuteration of water varies widely, depending on location and 
evolutionary time. At 1 AU it ranges from 0.001 for 0.1 - 0.2 Myrs, to 0.1 for 
0.3 Myrs and to 0.03 for times later than 0.4 Myrs. A minimum HDO/H 2 O ratio 
near 0.001 extends from the radial region around 0.1 AU at early collapse times 


to nearly 9 AU at later times. The predicted HDO/H 2 O ratios of Wakelam et 


al. 


al. 


(2014) and Tornow et al. (2014b) are slightly lower than those of Aikawa et 


(2012) but larger than the ratios observed by [Persson et ah (2014). An over¬ 
estimation of the water deuteration may be caused by a low collapse temperature 
which would suggest perfect cooling. Including the ortho-to-para ratio of H 2 can 
also lead to lower HDO/H 2 O ratios ( |Taquet et al.|[20T3 2014). In addition, water 
formed during the quasi-stationary core stage must be destroyed by the UV radia¬ 
tion of the interstellar radiation field, nearby stellar sources or X-rays (producing 
secondary UV photons in a similar way as cosmic rays) and UV photons from 
the protostar ( France et al.|20f4 ). As a result, photodissociation of deuterium-rich 
water molecules (HDO + hi/ —> D -|- OH, or H -|- OD), together with 1-D radial or 
vertical mixing enable a subsequent reformation of water by hot neutral reactions 
in the inner disk plane (Furuya et al.|[2M3). 
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R [AU] 


l.Oe+00 


l.Oe-01 


l,0e-02 


l.Oe-03 


l.Oe-04 


Fig. 15 HDO/H 2 O ratio in the gas phase on Lagrangian paths. At 1 AU the HDO/H 2 O ratio 
varies between 10“^ a nd 0.1, while in the du st phase it is much smoother, and only varies 
between 0.03 and 0.05 ||Tornow et al.|2014a|b[|. 


3.2 Fractionation of carbon and nitrogen in the solar nebula 


Isotopic enrichment of in cometary organics (the nitriles CN and HCN) 
and ammonia is commonly observed (Jehin et al. 2009[ Bockelee-Morvan et ah, 
this issues). Chemical models have shown that fractionation in low-temperature 
ion-molecule reactions can produce large enrichments in interstellar clouds 
( Rodgers and Charnley|[2008 Wirstrom et al. 20121. However, a recent study by 
Roueff et al. (20151 has demonstrated that several key reactions may in fact pos¬ 


sess energy barriers, leading to doubts as the efficacy of ion-molecule processes in 
N fractionation. Of particular interest for comets is that fractionation in the 
outer nebula should be similarly impaired by low values of the cosmic ray ioniza¬ 
tion rate, (. Thus, as found by Cleeves et al. (2014a I for water deuteration, the 
ratios commonly measured in both Oort Cloud and Jupiter Family 


low 


comets is strongly suggestive of an interstellar origin. 

The ^^C/^^C ratios measured thus far in simple cometary molecules (CO, 


HCN, CN) are solar (Jehin et al.|2009| and Bockelee-Morvan (these proceedings). 


tal ratio can increase to order unity ( 

Oberg et al.||2011 

. Because water, OH and 

O 2 recycle carbon nuclei back into CO ( 

Langer et al. 

1984), ion-molecule frac- 


tionation in regions with C/0~ 1 will act instead to channel carbon nuclei into 
simple organic molecules. Thus, as the so le carbon fractionat ion process involves 
the ^^C^ exchange reaction with ^^CO (Langer et al. 


19841 


CO condensation 

may in fact be the reason for the almost uniformity of cometary ^^C/^^C ratios 
measured to date. On the other hand, the condensed CO is predicted to produce 
distinctive ^^C/^^C fractionation patterns in the organic molecules formed in the 
grain-surface reaction scheme of Figure]^ ([Charnley et al. 2004). These involve 













































26 


K. Willacy et al. 


molecules whose isotopologues have not yet been detected in comets (e.g., 
H 2 CO, CH 3 OH) but for which the predicted fractionation has been confirmed by 
observations of protostellar envelopes (Wirstrom et al. 2011 1 . In this case, once 
again, ion-molecule pathways to fractionation in the nebula are predicated on 
the presence of significant flux of cosmic rays. 

If ion-molecule reactions indeed play a greatly reduced role in the cooler parts 
of the disk then the major sources of nebular fractionation will involve UV pho¬ 
tochemistry. Isotope selective photodissociation of CO has been proposed as an 


expla nation of the oxygen isotope ratios measured in meteorites (Lyons and Young 
2005). Isotopic enrichment of and D though isotope-selective photodisso¬ 


ciation of isotopologues of CO, N 2 and HD has also been considered (Woods and 
Willacy|2009 Heays et al.|2014 Cleeves et al.|2014a). However, these processes will 


only occur in thin upper layers of the disk and so will be negligible at the midplane. 
Thus, to be viable, these processes requre efficient downward vertical mixing in 
the disk. On the other hand, ISP occurring in the over layers of the presolar core, 
where the interstellar UV field is attentuated, could in principle produce isotopic 
anomalies in the oxygen and carbon isotopes (Yurimoto and Kuramoto|2004 Lee 


et al.|[2008 Charnley and RodgerspOOQ I. This environment could also ultimately 


produce significant fractionation in nitrogen-bearing species but this has not yet 
been demonstrated quantitatively. 


3.3 Solar nebula chemistry and the origin of oxygen isotopic anomalies in the 
solar system 


Oxygen is a very common element that is a major constituent of many minerals 
and rocks. It exists in three stable isotopes: ^®0, and The fractionation 
between these isotopes is not consistent across all solar system samples, as might 
be expected if the oxygen isotopes were homogeneously distributed across the solar 
nebula, and if only mass-dependent fractionation processes were active. Figure 
is a three-isotope plot of the oxygen isotopic compositions of primitive solar sys¬ 
tem materials such as chondrites, chondrules, shows that their relative isotopic 


compositions are different compared to those of terrestrial samples ( McKeegan et 
aL]|2011 Yurimoto et al.|[2007 Clayto n||1993[) The terrestrial line is characteristic 


of mass-dependent fractionation, where is twice as fractionated as relative 
to ^® 0 , whereas the steeper, slope -1 line, of the primitive materials suggests it is 
the result of mass-independent fractionation. 

There are three possible mechanisms that could achieve this mass-independent 


effect (Yurimoto et al. 2007). Firstly the isotope anomalies could represent an 
inherited heterogeneity in solar system materials, where the solids are rich in 
and the gas in the rarer isotopes. This could occur either from nucleosynthesis in 
other stars ( Clayton et al.|[l973 l or from nuclear reactions triggered by energetic 
particles produced by the Sun or by galactic cosmic rays ( Lee|1978 |. These nuclear 
processes could produce fluctuations in the abundance. The lack of observed 
isotop ic anomalies in other elements suggests this mechanism is unlikely ([ Cluj^d^ 


2004 


1993|). Additionally, presolar grains are rare in meteorites (|Nagashima et al. 
Nittler]|2003| ). 


A second suggestion is chemical fractionation within the solar nebula. Dif¬ 
ferences in the reaction rates of symmetric molecules, e.g. compared to 
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Cometary water here? 
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Fig. 16 A three-isotope plot for oxygen in the solar system. The Earth and Moon lie along a 
fractionation line with slope 0.52. This line is consistent with fractionation of a single homoge¬ 
nous source, by processes that depend on the mass of the isotopes, so that is nearly twice 
as fractionated as Data from meteorites and the Genesis mission falls along a slope-1 line 
indicating they contain material that has undergone some kind of mass-independent fraction¬ 
ation. The Genesis results suggest that much of the rocky material in the inner solar system 
was enriched in and relative to before the accretion of the first planetesimals. 
Hence rocky planets are not made of the average nebular material that formed the Sun. From 
[McKeegan et aL]||2011[| (reprinted with permission). 


asymmetric ones, e.g. can lead to the observed isotopic distributions, since 

the most abundant isotope is almost entirely contained in the symmetric molecules, 


while the rarer isotopes are in the asymmetric ones (Heidenreich and Thiemens 


1983). However, the exact reactions by which this effect might be achieved are 
uncertain (Yurimoto et al.||2007) 


The final possibility is isotope-selective photodissociation of CO, either in the 


parent cloud ( 

Yurimoto and Kuramoto 

2004 

) or at the surface of the protosolar 

disk 1 

Young|2007 Lyons and Young|200f: 

) could also result in enhancements of the 


rare isotopes of oxygen. CO dissociates when irradiated by far-ultraviolet photons 
with wavelengths between 91.2 nm and 110 nm. The dissociation process occurs 
in two steps, with CO first entering a bound excited state before dissociating. As 
a result the CO absorption spectrum consists of a large number of narrow lines at 
specific wavelengths determined by the vibrational and rotational levels involved. 
The change in mass between C^®0 and its isotopologues C^^O and C^®0 slightly 
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shifts the lines in the absorption spectra, and consequently the spectra for the 
different isotopologues do not overlap. Hence when CO is dissociated, either in 
the protostellar disk or in the molecular cloud core, the dif ferent isotopologues 
do not shield each other (see, for example, Visser et al.||2009 ). will become 

optically thick closest to the surface, in a region where and can still 

dissociate. When coupled with turbulent mixing that can transport the and 
towards the colder, shielded midplane of the disk this can result in the formation 
of water molecules that are enhanced in the rarer isotopes and which could later 
be incorporated into CAIs and other meteoritic material. 

The in situ measurements that Rosetta will make will directly probe the iso¬ 
topic composition of the volatile component of a comet. This will test the viability 
of the potential causes of the oxygen isotopic anomalies and provide new clues as 
to the origin of the cometary material and the conditions under which it formed. 


4 The origin of the nitrogen deficiency in comets 


Nitrogen is one of the most abundant elements in the PSN l Asplund et al.|[2009 l. 
Its abundance was measured in the giant planets atmospheres in which its stable 
thermodynamic form is NH 3 . The nitrogen abundance is found moderately to 


substantially enriched, compared to its solar value, in Saturn and Jupiter (Mousis 
et al. 20141, but is loosely constrained in Uranus and Neptune (de Pater et al.|1989 


de Pater and Richmond||1989 Ali-Dib et al.|[2M4 1. A puzzling revelation was the 


apparent severe depletion of nitrogen in comets. Mass spectrometry measurements 
of IP/Halley’s dust composition allowed to find that the nitrogen-to-oxygen ratio 
is depleted by a factor ~3 with respect to the solar value iJessberger and Kissel 


1991). Subsequent analysis showed that the ice part of the comet is depleted in 


nitrogen by a factor up to 75 (Wyckoff et al. 19911. This apparent depletion appears 


controversial because the N 2 molecule, thought to be the main nitrogen-bearing 
volatile in the whole PSN (Lewis and Prinn 1980 Mousis et al. 2002), has a 


molecular mass close to that of the very abundant CO molecule, implying that the 
mass spectrometer aboard the Giotto spacecraft did not have enough resolution 
to separate the two molecules ( Eberhardt et al^|1987 |. 

More recently, Cochran et al. (2000) used the 2DCoude spectrograph on the 


2.7-m Harlan J. Smith telescope of McDonald Observatory to observe the N^/CO 2 
in comets 22P/1995 SI (deVico) and C/1995 01 (Hale-Bopp), and found values 
around ~ 1 x 10“^. These values hint at a severe depletion in N 2 since the N^/CO^ 
ratio is usually equal to N 2 /CO. This is surprising specially that N 2 and CO have 
very close condensation temperatures ( Fray and Schmitt||2009 ). Moreover, efforts 
were made to find N 2 using the Far Ultraviolet Spectroscopic Explorer, but with 
no success (Bockelee-Morvan et al. ]2004l |. All these observations, and many more, 
motivated theorists to try to explain this depletion. 

Iro et al(|(2003), the authors assumed that volatiles are accreted by comets 


In 


as clathrates (jLunine and Stevenson||1985 ), and explain the nitrogen depletion by 


the low efficiency of N 2 trapping in water ices. One of this work’s conclusions is 
that N 2 can be trapped only if water ice is present in abundance higher than 2.8 
X solar value. This is due to the relative ease by which CO can be trapped in 
comparison with N 2 , so this last can be clathrated only if there is enough water 
left after all the present CO was captured. Since N 2 is trapped in crystalline 
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water ice at temperatures lower than CO ( Mousis et al.|[20T0 ), their work implied 
that comets might have formed at temperatures higher than 45 K, allowing the 
clathration of CO but not N 2 . This temperature lower limit is one of shortcomings 
of this model since in many protoplanetary disks models, temperature is thought 



ratio in clathrates since the potential parameters describing the guest-clathrate 
interaction in their statistical model correspond to guest-guest interactions (case 
of pure solutions) and not to the guest-water interactions usually used in such 


models (Mousis et al. 

2012 bI. Also, the model is unable to explain the thick N 2 

ice cover on Triton and Pluto, both bodies thought to be formed in the comets 

region ( 

Lellouch et al.||2011 

1 . 


Motivated by the shortcomings of Iro et al. (2003), Mousis et al. (2012b) in¬ 


troduced another statistical model based on an appropriate set of potential pa¬ 
rameters and that seems to address these issues. In this model, volatiles are also 
accreted as clathrate-hydrates, except N 2 which, due to its poor trapping propen¬ 
sity, is assumed to condense at ~20 K in the PSN and subsequently trapped in 
pure cristalline form in comets. Later radiogenic heating by the decay of short¬ 
lived nuclides will heat up these bodies to temperatures enough to liberate the 
trapped N 2 , but not CO. In Triton and Pluto, the gravity is sufficiently important 
to prevent the gaseous N 2 from escaping, and let it condense on the surface and 
form the observed cover. In the case of comets however, the gravity is negligible 
and the gaseous N 2 will escape, leaving comets depleted in nitrogen. A shortcom¬ 


ing of this model is that no radiogenic particles were found by stardust (Flynn 


2008). 


A recent model, applied to the description of Uranus and Neptune formations, 


also provides an alternative solution. Ali-Dib et al. (2014) proposed that the dif¬ 


fusive redistribution of vapor will remove the gaseous N 2 from the region below 
its iceline, possibly explaining the probable nitrogen depletion in Uranus and & 
Neptune if they formed in this region. One can invoke a similar argument for the 
nitrogen depletion in comets. In this case, the CO-rich comets will need to form 
between the CO and N 2 icelines, in order to accrete CO in solid form (since its 
vapor will also be depleted below its iceline) but not N 2 . The relatively small dis¬ 
tance between the two icelines (~4-5 AU) is however a caveat for this model since 
it imposes the formation of all comets in this restricted area of the outer PSN. 


Finally, all these models took the N 2 depletion at face value. As discussed 
above, the N 2 molecule’s mass is very close to CO, and the two molecules can be 
distinguished only using an instrument with a very high spectral resolution, such 
as ROSINA on board the Rosetta spacecraft ( Gulkis and Alexander|[2008 ). If this 
(or another) instrument is able to detect important amounts of N 2 in a comet, this 
might resolve at least part of the problem. Since important chemical differences are 


known to exist between Oort Cloud comets and Jupiter Family comets (Mumma 
and Charnley|20TT ), positive detections or dehnitive exclusion of N 2 in both types 
of populations are needed before drawing hnal conclusions. 
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5 Spin temperature as a cosmogonic indicator in comets 


The volatile composition, D/H ratio in molecular volatiles and the spin tempera¬ 
ture, as given by the ortho-to-para ratio in water or other molecules, are considered 
to be possible cosmogonic indicators for comets and as such may be interpreted 
as providing information about the formation history of the molecules in question. 
In this section, we discuss spin temperature as a possible cosmogonic indicator for 
comets. 

Molecules that contain identical hydrogen nuclei (like H 2 O and CH 4 ) display 
nuclear spin isomers. Water, for example, has two spin species referred to as ortho- 
H 2 O if the nuclear spins of the hydrogen atoms are parallel and para-H 2 if they are 
anti-parallel. The lowest energy level of para-H 20 lies 23.8 cm“^ (~34 K) below 
the lowest ortho level, so when water is formed at temperatures below about 50 K, 
the formation of para states is preferred and the ortho-para ratio (OPR) is less 
than the statistical equilibrium value of 3. Likewise for methane, statistical equi¬ 
librium corresponds to abundance ratios A:E:F of 5:2:9 and statistical equilibrium 
is reached for temperatures above about 60 K. 

OPRs have been routinely measured in water in comets since the apparition of 
comet IP/Halley | Mumma et al.|1987 ). This is currently done by fitting synthetic 
fluorescent emission models independently for each spin state as shown in Fig¬ 
ures and Corresponding spin temperatures are determined by placing the 
OPR on a theoretical curve such as that for water shown in Figure[^ OPR is also 
commonly inferred for ammonia by measuring the OPR for NH 2 , thought to be 
formed primarily via photodissociation of ammonia ( Kawakita et al.||2004 Shin- 
naka et al. ]WT| ). More recently, spin temperature has been measured in cometary 


methane for comets C/2007 N3 Lulin ( Gibb et al.||2012 |, C/2007 W1 (Villanueva 


let al.||2011[ |, C/2004 Q2 (|Bonev et al.|[2009 ' |Kawakita and Kobayashl] 
C/2001 Q4 ( Kawakita^et al.||2005 ). 


20091, and 


The meaning of spin temperatures has been the subject of much debate. For 
an isolated molecule of water or methane in the gas phase, for example, nuclear 
spin conversion is strongly forbidden. For this reason, spin temperature has been 
considered a cosmogonic indicator in comets that can provide information on the 


formation temperature (Crovisier 1984 Mumma et al. 19871. Observational evi¬ 


dence supports the suggestion that OPRs do not change once volatiles are released 
into the gas phase. For example, recent results mapping the spin and rotational 
temperatures in cometary comae are consistent with a scenario in which the OPR 
is unchanged by coma processes ( Bonev et al.|2068 |. Also, the OPR for NH 2 does 
not correlate with the heliocentric distance (and hence rotational temperature) of 
the comet, nor does it correlate significantly with abundances of other volatiles 
( Shinnaka et al.||2011 ). 

However, it has also been suggested that the spin temperature of water depends 
on the evaporation history of the ice ( Buntkowsky et al.|2008 |. Hama et al. (20111 
found experimentally that the OPR of water desorbed from ice after being vapor 
deposited at 8 K or formed via irradiation of a CH 4 /O 2 mixture at 8 K was consis¬ 
tent with statistical equilibrium, suggesting that OPR may equilibrate during the 
thermal desorption process or that energetically produced water does not undergo 
nuclear spin conversion (to a lower OPR) in water ice. Both formation processes for 
water in the Hama et al. ( |2011 1 experiment were energetic, and the authors note 
that the spin temperature of water molecules produced by non-energetic processes 
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Fig. 17 Spectrum of comet C/2007 N3 (Lulin) on 31 Jan 2009 with the sum of the telluric 
and best-fit synthetic fluorescent emission model (red) overplotted. The best fit water ortho 
(orange) and para (purple) models are shown. Also shown is OH prompt emission (olive). To 
illustrate the quality of the fit, the residual (scaled by a factor of two) with the best-fit models 
subtracted is shown. The 1-sigma error envelope is overplotted in green. The best fit is for 
statistical equilibrium (T^pj^ > 60 K). 
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Fig. 18 Residual spectrum of comet C/2007 N3 (Lulin) on 31 Jan 2009 (C) with the sum 
of the best-fit synthetic fluorescent emission model (red) overplotted (A) and offset vertically 
for clarity. The OH prompt emission is shown in the inset (B). To illustrate the quality of 
the fit, the residual (scaled by a factor of two) with the best-fit models subtracted is shown 
in (D) with the 1-sigma error envelope overplotted in green. Beneath the residuals, individual 
fluorescent emission models are shown for the three spin components of CH 4 (E) and CH 3 D 
(F) expanded vertically for clari ty. The best fit is for statistical equilibrium (Tspj^ >60 K). 
Adapted from [^bb et al^ | |2012| |. 
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H^O Spin Temperature (K) 

Fig. 19 OPR vs Spin Temperature in water with measurements for individual comets. In¬ 
dividual OPR measurements are from Mumma et al. (1988:1a, lb, 2); Crovisier et al. (1997: 
3, 1999: 4a, b); Dello Russo et al. (2005: 5, 6, 7a, 7b, 2007: 11, 12, 2008: 16); Kawakita et 
al. (2006: 8); Bonev et al. (2009: 9, 2008: 10); Woodward et al. (2007: 13); Villanueva et al. 
(2011: 14); Paganini et al. (2014: 15, 2015: 17, 2012: 18); DiSanti et al. (2013: 19); Radeva et 
al. (2010: 20); Gibb et al. (2012: 21). 


must also be studied before definitive conclusions can be drawn, especially since 
cometary volatiles may have formed by low-energy processes. Nonetheless, such 
studies raise the possibility that cometary OPRs may not, in fact, be cosmogonic. 


Can the measurements of OPR in comets shine any light on this question? 
We show the measurements of cometary OPR in water to date in Figure |19[ 
Water was chosen due to the greater number of OPR measurements than for other 
molecules in comets. Early measurements suggested a clustering near T^pin ~ 30 K, 
leading many authors to conclude a formation temperature of ~30 K for a large 
population of comets. Similar conclusions have been drawn from studies based on 
CH 4 ( Kawakita et al.|2006 ) and NH 2 i Shinnaka et al.|2011[ |, particularly given the 
lack of correlation between spin temperature and other properties of the observed 
comets ( Shinnakaet al.|2011 1 . However, it must be noted that among the 20 comets 
measured to date, nearly all of them are within 1 - or 2 -sigma of the statistical 
equilibrium value, excepting the measurements reported for IP/Halley, one of the 
two measurements for C/2001 A2, and C/1995 01 (Hale-Bopp). Likewise, half 
of all measurements are at or above the statistical equilibrium value. It is also 
notable that no measurements exist below a spin temperature of 20 K, and this 
observation holds for methane and ammonia spin temperatures as well. 


Possible interpretations of spin temperatures measured in comets are therefore 

( 1 ) that it is a cosmogonic indicator that gives the formation temperature of the 
molecule in question and therefore perhaps the formation conditions of the comet, 

( 2 ) that it is a measurement that, while unchanged by coma processes, is set by 
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processes occurring in the solid ice during or prior to sublimation. Current data 
are not sufficient to draw firm conclusions at this point. Even if further observa¬ 
tions result in most comets consistent with statistical equilibrium, we may still 
be sampling cosmogonic material since it is known that comets contain material 
processed at high temperatures, such as crystalline silicates, which are interpreted 


to originate from turbulent radial mixing in the solar nebula (Bockelee-Morvan 
et al.|[2002). Mixing of volatile material that formed at high temperatures in the 


inner disk to the comet forming region would have resulted in a higher OPR than 
for uncontaminated material that formed in the cold midplane. Clearly a great 
deal more work needs to be done to understand the significance of molecular spin 
states in cometary material. 


6 On the origin of neutral sodium in comet tails 


The presence of sodium D line emission has been confirmed in a large number 
of comets close to perihelion since it was first reported more than a century ago 
( Bredichin[1882 ). Observations of comet C/199501 Hale-Bopp during the spring of 
1997 led to the discovery of a new tail connected with the sodium D line emission. 
This neutral sodium gas tail is entirely different from the previously known ion and 
dust tails, and its associated source is unclear. Up to now, only physical reasons 
have been advanced to rationalize its origin. It has been proposed that this third 
type of tail is shaped by radiation pressure due to resonance fluorescence of sodium 


atoms (Cremonese et al. 1997). Further possibilities were considered and then 


rejected, as photo-sputtering and/or ion sputtering of nonvolatile dust grains (Ip 
and Jord^|1998 1 . In the same way, collisions between the cometary dust and very 
small grains were also questioned ( |Ip and Jorda|[l9M I. The scenario presented in 
this review is completely different since entirely based upon chemical grounds. It is 
shown that the Na+ ions, washed out of the refractory material during the hyration 
phase, loose their positive charge to evolve progressively into neutral species with 
their migration towards the surface of cometary ices. 


Ellinger et al. (2015) follow the chemical path of sodium, starting from the 


incorporation into the ice until the final transformation into a neutral atom when 


released from the sublimating crystalline (see Section 2.5) cometary ice. The mod¬ 
eling of bulk, surface structures and reactivity of ices has been essentially developed 
for environmental applications ( Casassa and Pisani]200^ Casassa et al.|2d05 Ca- 
latayud et al. 2003) Levering on these results, high-level numerical simulations 


based on first principle periodic density functional theory (DFT) were used to de¬ 
scribe the solid structure of the ice in the form of apolar hexagonal ice (Bussolin 
et al.|[T998 l composed of bi-layers of water molecules. 

A code, specifically designed for the study of periodic systems, namely, the 


Vienna ab initio simulation package (VASP) (Kresse and Hafner 1994a|b ) was 
used to carry out all the calculations. The generalized gradient approximation 
(GGA) was employed in the form of the Perdew and Wang exchange-correlation 
functional ( Perdew et al.|[T992 | coupled to the Grimme correction (Grimme et al. 
2010) to take care of the long-range van der Waals interactions. The evolution of 


the ionic character of the sodium atom as it reaches the ice surface was obtained by 
a topological analysis (Silvi and Savin|1994) performed by means of the TOPMOD 
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package | Llusar et al. 1999|| and its recent extension to periodic systems (|Kozlowski 
and PilrnepOll 


A first calculation on isolated Na(H20) showed that the interaction of Na with 
one H 2 O is similar to that between two H 2 O molecules in the H 2 O dimer (« 0.3 
eV), which gives a hint in view of a possible replacement of a molecule of water by 
a sodium atom in organized structures. In liquid water, microsolvatation of sodium 
atoms has been studied on small clusters of Na(H 20 )n (n=l-6) to mimic the water 
environment. Both the structures in which Na is positioned on the surface and the 
encapsulated structure in which Na is surrounded by a solvation shell are minima 


on the potential surface. They are very close in energy (Hashimoto and Morokuma 


1994) 



Fig. 20 Typical structures of Na-H20ice arrangements. Top: Na in the cavity inside hexagonal 
ice. Middle: Na replacing H 2 O in the crystal. Bottom: Na adsorbed on the surface of the ice. 
Color code Yellow (Na); red (O) ; white (H). 


In the ice, the periodic structure of the crystal imposes strong geometrical 
constraints. The most stable structure, not found in the liquid, is the one where 
Na takes the place of a water molecule in the lattice (middle structure in Figure [M| 
when the ice is formed (Z\E « 1 eV). Contrary to microsolvatation, the stability 
of Na encapsulated in the bulk of the crystal {Z\E « 0.5 eV), shown in the top 
structure in Figure 201 differs from that of the adsorbed Na according to the 
position on the surface. For example, the adsorption energies vary from (AE « 
0.6 to ^ 0.5 and si 0.1 eV) when the adsorption takes place over the center of a 
surface hexagon, over a surface oxygen (bottom structure in Figure 20) and over 
an OH dangling bond, respectively. 

The most relevant point to this study is that the charge on the sodium 
atom depends strongly on the position with respect to the ice surface, as does the 
adsorption energy. When constrained to stay in the cavity inside the hexagonal 
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lattice, Na remains largely ionized {qNa ~ +0.9) which is close to the value of 
« +0.8 when Na takes the place of an H 2 O molecule in the bulk. When Na is 
adsorbed on top of an oxygen of the surface layer, qNa~ +0.2, a value similar 
to that found for the isolated Na(H20) complex. Then, as the ice is progressively 
eroded by the radiation field, Na eventually reaches the surface where it can desorb 
as a neutral atom when the upper layer of the ice vaporizes. 


7 Summary 


The observed compositions of comets are the end product of a chemical history 
that began in the cold parent molecular cloud core. The change in physical condi¬ 
tions experienced by the core material as the core collapsed and the protosun and 
protosolar disk formed and evolved led to chemical changes that are reflected in 
the cometary volatiles that we observe today. Each stage of evolution left its im¬ 
print in the mixing ratios of molecules observed in comets and disentangling their 
effects requires interaction between comet and disk researchers. The combination 
of the two fields of study provide a powerful tool in understanding the drivers 
that determine comet composition and conversely what the variations in comet 
composition can tell us about the relative formation conditions of these comets. 

Protostellar disks provide complementary information to the comet observa¬ 
tions. Their chemistry is complex, covering ice chemistry in the cold, shielded 
midplane, warm gas phase chemistry and hot, highly ionized chemistry. Comets 
form from the ices in the midplane, yet dynamical processes ensure that this region 
does not evolve in isolation but rather is the end product of material that has been 
processed in a variety of different environments and then brought together in the 
cometary ices. Establishing the origin of cometary ices therefore requires not only 
studying the midplane but also an understanding the chemical and physical evo¬ 
lution of the disk as a whole. By studying disks we can observe the chemistry and 
physics that shaped the solar system as they shape other protoplanetary systems. 

Molecular deuteration is a link that might tie together the disparate objects 
of the solar system and their formation history and give us the picture of the 
whole solar system formation history. However, at present this link has still weak 
connections that need to be reinforced. The advent of the powerful new facilities, as 
for example ALMA and IRAM-NOEMA, will certainly allow this field to progress 
rapidly in the next few years. 

Rosetta will provide detailed ground-truths for the organic and isotopic makeup 
of cometary volatiles. The new insights obtained into the chemical processes which 
likely produced these molecules - interstellar or nebular - will allow us to constrain 
models of the formation and evolution of the early solar system. Its instruments 
will characterize the composition of the comet 67P/Churyumov-Gerasimenko and 
distinguish between the molecules native to the comet nucleus (parent species) 
and those generated by chemistry or by fragmentation in the comae (daughter 
species). The observations it will make will provide new insights into the formation 
conditions of the comet. They will also elucidate the degree to which interstellar 
material can survive its journey into the solar nebula and from there into comets. 

Several measurements will characterize the its thermal history, e.g., measure¬ 
ments of molecular spin temperatures ( Gulkis and AlexanderpOOS I and deutera¬ 
tion and of the abundance of argon (Hassig et al.||20T4|. The relative abundances 
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of elements can tell us something about the dominant molecular form of each el¬ 
ement at the time of formation - if an element is depleted relative to others then 
the major carrier of that element must have been too volatile to be condensed at 
the time of formation, or must have desorbed in the intervening time. In particular 
Rosetta measurements may help to elucidate the mystery of the lack of nitrogen 
in comets ( Altwegg et al.||2Q15 |. 

A key measurement will be the ratio of oxygen isotopes in water. Meteoritic 
water has a very different ratio of the oxygen isotopes compared to oxygen isotopes 
in the Sun. Rosetta provides our best opportunity to directly detect the isotope 
ratio from the early solar system, and in so doing will distinguish between the 
various theoretical models of how the oxygen isotope fractionation occurred. 

In summary, Rosetta’s observations will provide new insights into the chemical 
processes that resulted in the formation of cometary volatiles and consequently into 
the conditions that were present in the solar nebula at their time of formation. 
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